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Abstract 

Fermions coupled to Yang-Mills matrix models are studied from the point 
of view of emergent gravity. We show that the simple matrix model action 
provides an appropriate coupling for fermions to gravity, albeit with a non- 
standard spin connection. Integrating out the fermions in a nontrivial geomet- 
rical background induces indeed the Einstein-Hilbert action and a dilaton-like 
term, at least for on-shell geometries. This explains and precisely reproduces 
the UV/IR mixing for fermions in noncommutative gauge theory, extending re- 
cent results for scalar fields. It also explains why some UV/IR mixing remains 
even in supersymmetric models, except in the N = 4 case. 
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1 Introduction 

Recently it was understood that gravity can arise effectively from non- commutative gauge 
theory [1-3]. This mechanism is realized in certain matrix models of Yang-Mills type, 
which are known to describe gauge theory on non- commutative (NC) spaces. These 
models can be interpreted more naturally as describing dynamical NC spaces, suggesting 
an intrinsic realization of some sort of gravity theory. The mechanism of this "emergent 
gravity" is quite simple and intrinsically non-commutative [1]: On a general background 
corresponding to a generic NC space, all fields in the model couple universally (up to 



possibly density factors) to an effective metric or frame associated with the background. 
The Einstein-Hilbert action is induced upon quantization. Such a generic NC background 
can be equivalently interpreted in terms of a U(l) gauge field on a fixed Moyal-Weyl space 
M.q. This then provides an understanding of UV/IR mixing in the one- loop effective action 
of NC gauge theory, in terms of an induced gravity action. 

In the present paper, we include fermions in these matrix models, and study their cou- 
pling to the emergent gravity. In particular, we extend the explanation of UV/IR mixing 
in terms of induced gravity [4] to the case of fermions. The matrix-model framework 
strongly suggests a simple action for fermions, as realized e.g. in the string-theoretical 
IKKT model. We study this fermionic action in detail on a generic NC background, fol- 
lowing the geometrical point of view of [1]. This leads to an effective semi-classical action 
for a fermion on a background with metric G a b(y). The action is similar to the standard 
action for fermions on a curved background, except that the spin-connection appears to be 
missing. It is nevertheless well-defined, due to the existence of an (unobservable) preferred 
frame, which is an intrinsic part of emergent NC gravity. 

The main result of this paper is that in spite of this unusual feature, the simple 
fermionic action under consideration is nevertheless reasonable on general backgrounds, 
and should be physically viable. In the point particle limit, the fermions will follow the 
same trajectories as for the conventional coupling to gravity, albeit with a different rota- 
tion of the spin along the trajectory. Furthermore, we determine the effective gravitational 
action obtained by integrating out the fermions in our framework. This leads indeed to 
the expected Einstein-Hilbert action for the effective metric, with an extra term for a 
density factor, and another term which vanishes on-shell. We conclude that the fermionic 
action considered here is the "correct" one for emergent gravity, suitable for a physically 
realistic theory of gravity. 

For a consistent quantization of emergent gravity, a supersymmetric extension of the 
model appears to be necessary. We point out that the induced gravitational action will 
be finite if and only if the model enjoys N = 4 supersymmetry at the Planck scale. This 
explains in particular the persistence of UV/IR mixing in NC gauge theory in models with 
less supersymmetry. Such a supersymmetric extension is realized by the IKKT model, 
which is known to admit noncommutative backgrounds; see e.g. [5-8] and references 
therein for related work including (indirect) evidence for gravity on such backgrounds. 

This paper is organized as follows. In Section 2, we review the basic aspects of emer- 
gent gravity, and write down the actions under consideration. In Section 3, we review 
the mechanism of induced gravity, and set up the form of the appropriate Seeley-de Witt 
coefficients corresponding to the fermionic action. They are the ingredients which pro- 
vide the induced effective gravity action, and they are modified in comparison with the 
commutative case. These coefficients are computed in Section 4, which is the core of this 
paper. We also provide an exact expression for the Ricci scalar in terms of the effective 
metric of emergent gravity; this should be useful for other considerations as well. In 
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Section 5, we consider the same problem from the point of view of NC gauge theory, and 
rewrite the obtained gravitational action in terms of U(l) gauge fields on Moyal-Weyl 
space. The one-loop effective action for this NC gauge theory is carefully computed in 
the appropriate IR limit, which provides a precise matching with the gravitational action 
obtained before. We conclude with a short discussion of supersymmetry and an outlook. 



2 Matrix models and effective geometry 

As a starting point of emergent NC gravity, we consider the semi-classical geometry arising 
from the matrix model with action 

S YM = -Tr[Y a ,Y b ][Y a ',Y b '}g aa/ g w . (1) 

Here Y a , a = 1,2,3,4 are hermitian matrices or operators acting on some Hilbert space 
Ti,, which constitute the dynamical objects of the model. The model also contains a 
constant background metric, which is 

g aa > = 5 aa > or g aa > = i] aa i (2) 

in the Euclidean resp. Minkowski case. The commutator of 2 matrices is denoted as 

[Y a , Y b ] = id ah , (3) 

defining 9 ab e L{H) as an antihermitian operator-valued matrix. We focus here on con- 
figurations where the Y a (which need not be solutions of the equation of motion) can be 
interpreted as quantization of coordinate functions y a on a Poisson manifold (Ai,8 ab (y)) 
with general Poisson structure 9 ab (y). This defines the geometrical background under 
consideration, and conversely any Poisson manifold provides (locally) after quantization 
a possible background Y a . More formally, this means that there is an isomorphism of 
vector spaces 

C(M) -> A C L(H) 

f(y) ^ f(Y) (4) 
i{f,g) [f,g} + 0(9") 

where C(M) denotes the space of functions on Ai, and A is the algebra generated by 
Y a , interpreted as quantized algebra of functions. Thus 6 ab reduces in a semi-classical 
limit to a classical Poisson structure 8 ab (y). Observe that the Y a correspond to preferred 
coordinate functions y a on A4, and g a b = 5^ resp. g a b = r] ab defines a flat "background" 
metric, which is constant in the preferred coordinates y a . Indices throughout this work 
will always refer to these preferred coordinates^, and only some (final) formulas will be 
manifestly covariant. 

1 this is the reason why many of the formulas in this paper are written in a non-covariant way. A 
covariant formulation is possible but will not be given here. 
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To simplify things, we restrict ourselves to the purely geometrical or "irreducible" case 
in this paper, i.e. we assume that A is in some sense dense in L{H). This means that any 
matrix in L{7i) can be considered as a function of Y a resp. y a . From the gauge theory 
point of view in section |5], it means that we restrict ourselves to the U(l) case; this is 
most interesting here since the UV/IR mixing is restricted to the trace-£7(l) sector. For 
the nonabelian case see [1]. 

Besides the preferred coordinates y a , the noncommutative background provides in the 
semi-classical limit a preferred frame 

given in terms of the antisymmetric Poisson tensor, which we assume to be non-degenerate. 
This formula is only valid for the preferred coordinates y a , and does not admit local 
Lorentz transformations. The effective metric arising from the matrix model turns out to 
be [1] 

G ab (y)=9 ac (y)6 bd (y)g cd , (6) 

which has the correct tensor structure and is associated to the above frame. G ab (y) is 
indeed the effective gravitational metric in the matrix model (up to a rescaling discussed 
below), because it enters the kinetic terms for matter fields through [Y a , ty] ~ i9 ab (y)-^ty; 
this will be seen explicitly below. Note that G ab (y) is not flat in general. Hence the back- 
ground M. naturally acquires a metric structure (A4, ab (y), G ab (y)), which is determined 
by the Poisson structure and the constant background metric g^. While this is rather ob- 
vious for scalar fields, it turns out that essentially the same G ab also couples to nonabelian 
gauge fields [1], and to fermions as we will show here. 

An infinitesimal version of the metric (^|) was observed already in [2]. The frame (|5]) 
was also pointed out in [3] , as well as a metric of type (|6]) for the self-dual case. There is 
also some overlap with the ideas in [9]. For further related work see e.g. [10-12]. 

Scalars. We first review the case of scalar fields i.e. hermitian matrices $ coupled to 
the matrix model (|l|). The only possibility to write down kinetic terms for matter fields 
is through commutators [Y a , $] ~ i8 ab (y)-^ E ^, and one is lead to the action 

S[$) = -(27r) 2 Tr W [F a ,$][r a ',$] 

~ Jd^ P {y)G«\y)^{y)^{y)- (7) 

Here and throughout this paper, ~ indicates the leading contribution in a semi-classical 
expansion in powers of 9 ab . This involves the symplectic measure 

p(y) = (dete ab (y))-^ = \G ab (y)\^ = e~° (= A% c (y) ) (8) 



3 



on (Ai, 6 ab (y)), which can be naturally interpreted as non-commutative scale A% c ; here 
\G ab {y)\ = \det G a b(y)\. After appropriate rescaling of G ab (y), this can be rewritten in 
covariant form 

S[*] = f d A yG a \y)d y «§d yb <$> (9) 

with the effective metric 

G ab = \G ab \ 1/4 G ab = p{y)G ab 
\G ah \ = 1 (10) 

which is unimodular in the preferred y a coordinates. Recall that these are characterized 
by the background metric being constant, g ab = 5 ab resp. g ab = r] ab . 



Fermions. Then the most obvious (perhaps the only reasonable) action for a spinor 
which can be written down in the matrix model framework^ is 

S=(2n) 2 W la [Y a ,^] ~ J d A yp{y)^ilaO a \y)d b ^> (11) 

ignoring possible nonabelian gauge fields here to simplify the notation. This is written for 
the case of Minkowski signature, the Euclidean version involves the obvious replacement 
^ — > W. This defines the (matrix) Dirac operator 

TP = la [Y a , ~ i la ab (y)d b *. (12) 

We can compare this with the standard covariant derivative for spinors (see e.g. [13]) 

^comm * = H° e£ (d^ + S fec U%) * (13) 

where 

= 4^( V m4) (14) 

is the spin connection, and 

i 

Safe = ^[7a,7fe] • (15) 

Comparing (|j"2"|) with flT3|), we observe again that in the geometry defined by (fjl), 

e${y) := ^(y)^ (16) 

plays the role of a preferred vielbein. However this must be used with great care, because 
the distinction between the coordinate index fi and the Lorentz index a is lost in the 
special "gauge" inherent in (|T6]). 

2 In particular, fermions should also be in the adjoint, otherwise they cannot acquire a kinetic term. 
This does not rule out its applicability in particle physics, see e.g. [14]. 
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One notices immediately that the spin connection appears to be missing in (111), which 



seems very strange at first. One of the main messages of this work is that in spite of 



this strange feature, the action (JTT ) is a good action for a fermion propagating in the 
geometry defined by G a t,. Recall that the spin connection determines how the spinors are 
rotated under parallel transport along a trajectory. However, the spin-connection u ab can 
always be eliminated (via parallel-transport resp. a suitable gauge choice) along an open 
trajectory. Then the conventional kinetic term ([H|) boils down to (|ll|). Therefore in the 
point-particle limit, the trajectory of a fermion with action flTT| ) will follow properly the 
geodesies of the metricQ G a b, albeit with a different rotation of the spin. Furthermore, we 
will show in detail that the induced gravitational action obtained by integrating out the 
fermion in (|TTD indeed induces the expected Einstein-Hilbert term J d A y R[G] A 2 at least 
for "on-shell geometries" , albeit with an unusual numerical coefficient and an extra term 
depending on o. All this shows that ( |TTD defines a reasonable action for fermions in the 
background defined by G^. 

In particular, the transport along a closed curve determines the holonomy in a gravi- 
tational field, and the missing spin connection in the above action strongly suggests that 
holonomies here will be different than in General Relativity. More generally, the gravita- 
tional spin rotation for a free-falling fermion might provide a nice signature for or against 
the emergent gravity framework. 



Equations of motion. So far we considered arbitrary background configurations Y a 
as long as they admit a geometric interpretation. The equations of motion derived from 
the action ([[]) are 

[Y a ,[Y a \Y b ]]g aa , = 0, (17) 
which in the semi-classical limit amount to 

9 ma d m 9 nb g ab = or G ac d c 6j = 0. (18) 

They select on-shell geometries among all possible backgrounds, such as the Moyal-Weyl 
quantum plane (|59|). In the present geometric form they amount to Ricci-flat spaces [1,2] 
at least in the linearized case. However since we are interested in the quantization here, we 
will need general off-shell configurations below. The equations (|18"D are only valid in the 
preferred y a coordinates as discussed above, and the same applies to most computations 
below. They can be cast in covariant form, which will not be done here. 

3 for massless particles, the geodesies of G a b coincide with those of G a b- Masses should be generated 
spontaneously, which is not considered here. 
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3 Quantization and induced gravity 

We are interested in the quantization of our matrix model coupled to fermions. In prin- 
ciple, the quantization is defined in terms of a ( "path" ) integral over all matrices Y a and 
In four dimensions, we can only perform perturbative computations for the "gauge 
sector" Y a , while the fermions can be integrated out formally in terms of a determinant. 
Let us focus here on the effective action I\ obtained by integrating out the fermionic 
fields, 



Later, we will consider an alternative interpretation as Laplacian of the fermionic fields on 
R~ coupled to an adjoint U{1) gauge field. In Feynman diagram language, (^) will then 
amount to the sum of all one-loop diagrams with arbitrary numbers of external A-lines. 

In [4], the analogous computation of the effective action T$ for scalar fields was carried 
out in two different ways. On the one hand r$ was evaluated from the geometric point of 
view via induced Einstein-Hilbert action. On the other hand the action ([?[) was regarded 
as NC gauge theory by means of covariant coordinates, where the effective action for 
the NC gauge fields was determined by integrating out the scalar fields. The obtained 
effective actions were shown to agree in the IR regime as expected. As a consequence 
UV/IR mixing of NC gauge theory can be interpreted as an effect of gravity. 

In this paper, we generalize the results of [4] to the case of fermions. The effective 
action is evaluated in two different ways: We first compute as induced gravity 
action, using its semi-classical geometrical form. This is then compared with the one-loop 
effective action of NC u(l) gauge theory. As expected, we find complete agreement in a 
suitable IR regime. 

Square of the Dirac operator and induced action. Starting from the action 




(19) 



which for non-interacting fermions is given by 



r* 




(20) 



S= (2tt) 2 Tr&"f a [Y a ,V] 



(21) 



we want to study 




(22) 
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at one loop, considering the Euclidean case for the sake of rigor. The square of the Dirac 
operator takes the following form 

^ = lalb [Y a ,[Y b ,n 

-iai b e ac d c {e bd d d m) 

= -G cd d c d d m - a d d<$>, (23) 

with 

a d = iai b e ma d m e dh = -2i T, ah e ac d c e bd + 9ab e ac d c e bd . (24) 

See Appendix C for a comparison of this term with the commutative case, Jp 2 defines the 
quadratic form 

Square = (2tt) 2 M» ^ ~ J d% P (y)¥ ^ 

dfiylG^Wp 2 *, (25) 



which is very similar to the scalar action. In terms of the metric G a b (0) with \G 



ab 



^square can be written in covariant form 



5 square = j d A y VJZ)^, (26) 

in terms of the rescaled squared Dirac operator 

ffiv = ~ (G cd d c d d V + e~ a a d d d ^ . (27) 
We now compute the effective action using 

1 . . 1 poo 1 

-Tr(log^) 2 -log^ 2 ) = — Tr / da- ( e~ a ^ - e~ a ^ ) 

2 2 J a V / 

= --Tr T-fe-^-e-^e-^, (28) 

where A 2 denotes the cutoflQ for \jp 2 -, regularizing the divergence for small a. Now we 
can apply the heat kernel expansion 

Tre-^ = V«^ / d 4 ya n (y,fl 2 ) (29) 
„>o J m V 7 



4 We write 2 A 2 instead of A 2 in (g8|) in order to be consistent with the cutoff for | A 2 for scalar fields 
used in [4] , which is implicit in T$ as given below. 
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where the Seeley-de Witt coefficients a n (y,^ 2 ) are given by [15] 

1 

16^ 2 



a (y) = 77-0 trl, 



f = -G m "^ m n n + n m n n -r^ fe ,j (30) 
n m = ld mn (e^a n + f n ) (31) 



2 

and tr denotes the trace over the spinorial matrices. The effective action is therefore 

r * = 16^ /A ^2tr(ll)A 4 + tr ^^lll + ^ A 2 + 0(logA)j , (32) 

where tr(l) =4 for a Dirac fermion. Everything is expressed in terms of the unimodular 
metric G a b, which can be written in terms of G ab using 

R[G] = p(y)^R[G] + 3A G a-^G ab d a ad b a 
A G a = -(G ab d a d b a-T c d c a) 

= pj^=(detG ab ) 1/4 

pa = Qcd^ = e -a T a _ Q ba _ (33) 



The actual computation of tr£ is given in Section |4.2| . Note the relative minus sign of 
the various terms in the effective action I\ compared with the induced action due to a 
scalar field [4], 

J d 4 y (-2A 4 - \r[G}1 2 + 0(log A)) . (34) 



16 7T 

hence 

r*+4T $ = — ?— [d 4 ytx£A 2 . (35) 
16 7l 2 J 

This expresses the cancellation of the induced actions due to fermions and bosons, apart 
from the £ term. For the standard coupling of Dirac fermions to gravity on commutative 
spaces, one has [16] 

^comm = R (36) 

which originates from an additional constant term —jR in p 2 omm (Lichnerowicz's for- 
mula). In our case, £ turns out to be somewhat modified due to the vanishing spin 
connection, nevertheless it contains the appropriate curvature scalar plus an additional 
term (1571). This will be discussed further in Section |572|. 
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4 Computation of the Seeley-DeWitt coefficient 



In this section we determine the second Seeley-de Witt coefficient for We first obtain 
an exact, compact result for the Ricci scalar R[G] expressed in terms of the Poisson tensor 
8 mn (y). This may be useful for other purposes as well. We then compute trS in terms of 
9 mn (y), which turns out to be closely related to R[G] as desired. 



4.1 Ricci scalar in terms of 9 mn 

The curvature is given as usual by 

Rabc = d b T d ac - d a T bc + T e ac T d b - T e bc T d a . (37) 

The Ricci scalar is then 

R = G ac R h abc = G ac (d b T b ac — d a T b bc + T e ac T b eb — r^X^) . (38) 
In terms of the metric and its derivatives R is given by 

R = (d b G bd ) G ca (d a G cd ) + G ab G cd d b d d G ac - G mn G™d p d q G mn 
- (d b G bd ) G mn (d d G mn ) - -G™ (d p G mn ) (d q G mn ) 

+^G np (d p G ac ) (d c G na ) - l -G™G mn (d p G mn ) G kl (d q G kl ) . (39) 

We use the explicit formula for the metric tensor 

G mn (y) = 6 ma (y)6 nb (y)g ab (40) 
to express R in terms of 9 (see Appendix A for details) 

R = -2 (d a 9 ap ) G bc (d c 6 b p) - (d a 9 ap ) (d b 6 bq ) g pq 

+2G m ^e m \d p d q 9~lg ab + l -G bc (d c B- v l ) G ad (d d 0^) g" 

+26 mn Gm p d q e m l - l -G mn G^ {d p e m \) {d q e-l) g °» 

+46 bc {d b 6 da ) (d d a) g ca + 1(1* (d p 6 mn ) (d q 6^ n ) 

_ Q m P {dp9m) _ 1 { Q nQCq) {dc0 n P) ^ ^ 

This equation holds in fact for any vielbein using the identification 

e» m g mb = e\ 0^9 mb = ~e b (42) 
9 



since we have not exploited any property of 9 mn yet. However, by making use 
Jacobi identity, 

Wc 1 + dcO-J + d h e-J = o 
d a e pq = - (d c e-^) (e mp e cq - e mq e cp ) , 

several terms appearing in the computation of tr£ and R[G] are equivalent^: 

G mk (d k e™) {d n e m \) = hs* (d p e mn ) (d g e^ n ) 

{9 mn d q e-l) 9 qa G pk (d k B£) = 2{d m 6 ma )G nk {d k d-l) 
(6 mn d q 9-l) 6« a (d p 6 pb ) g ab = 2{d m 6 ma ){d n d nb )g ab 

e ma d m dj nb g ab = G mn d m d n a + 6 ma (d m 6 nb )(d n a)g ab 

= \o mn G m d p d q e^ n + Ig« (d p e™) {d q e m l n ) 

Thus in our framework the Ricci scalar associated to G ab takes the form 

R [G] = e™G™d p d q e^ n + G™{d p 9™){d q e^ n ) 

+2 (d m 6 ma ) G nk (d k 6£) - G mn (d m a) (d n a) 

+\c mk {d k e- n l) G nl (d^l) g ab - ^-G mn G pq (9 P C) (d q Kl) J* 

-~ (d m e na ) (d n e mb ) 9ab . 

Evaluating also 

3A G a -^G mn d m ad n a = ~G mn {d m 6 m ) (dj£) - \& mn G pq d v d q Q^ n 

+~G mn d m ad n <j - 3 (d m 6 ma ) G nk (d k 9^) 

gives 

= e"* [ - ~9™G™d p d q 9^ n - lew (d p 9™) {d q 9-l) 
- (d m 6 ma ) G nk (d k e~l) + X -G mn (d m a) (d n a) 

5 By means of these relations one can also check that the action (|2^) is indeed hermitian. 
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+\o mk {d k e-l) G nl {d^l) g ab - l - G ™G m {d p e m \) {d q 9-l) g ab 

-\{d m e na ) {d n e mh ) g ab \ (48) 
Via partial integration, the number of independent terms can be further reduced: 

d 4 y e- a 9 ma d m d n 9 nb g ab = 0, 
d'y e~° (d m 6 na ) {dj mb ) g ab = -J d'y tl (6 mn G™d p d q 6^ n + G™ (d p 6 mn ) (d q B^ n )) , 

j d A y e-° {d r fin G^ {d k 6-J) = J d A y e~° (d m 9 na ) (d n 6 mb ) g ab . (49) 
This yields the following compact form for the Ricci scalar 

d*y R[G]A 2 = J d'y e- { l -G mk (d k 9^) G nl (d^) g ab - l -G™G™ {d p 6 m \) {d q &£) g ab 
-~ (d p 9n G" k (d k e^) + l -G mn (d m a) (d n a) }X 2 , (50) 



which is indeed what we need because det G ab = 1. However, one has to be careful when 
using partial integration. We regard here the cutoff A 2 as a constant effective cutoff for 
Ag, independent of y a . In Section |5|, we will relate A 2 with the effective cutoff A for A^, 
which is the Laplacian of NC gauge theory. This corresponds to a position-dependent 
A 2 as given in flBGf). In that case, one must either include additional terms from partial 
integration, or use ( fiSD as we will do. 



4.2 Evaluation of tr S 

We need to evaluate 

tr£ = -tiG ab (d a n b + tt a tt b - f ^fi r ) , (51) 



where 



= \ {G mnlalb 6 pa (d p 6 nb ) - G mn (d p G pn ) + d m a) (52) 
and d n = e~ a a n . For the computation of tr£ we use the relations given in Appendix A 



and the Jacobi identity (|45|) . Recalling 



tryY = 4g ab } 



tr 7 a 7 6 7 c 7 d = 4 (g ab g cd - g ac g bd + g ad g bc ) , 

(53) 
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we find for the individual parts of tr£: 
trG mn d m n n = 2e-°{G mn {d m G np )g ab 6« a {d q 6 pb ) 



+ (d m 6 qa ) {d q 6 mh ) g ab + e ma d m d n e nb 9ab 

-G mn (d m G np ) {d q G qp ) - d m d n G mn + G mn d m d n a 

2e-°{ (d rn e ma ) c nk {d h e-l) - e ma d m d n e nb 9ab 
A G ™{d p e mn ) (d q e m l n ) + \e mn G™d p d q e 



r 1 

= (54) 

tiG mn Q m Q n = e- ff {{g ab g cd -g ae g M + g ad ^)G nk e qa (d 9 e nb )e lc (d l e kd ) 
-2 (d q G« n ) G nk g cd 6 lc (dt9 kd ) + 26^ (d q 6 nb ) (d n a) g ab 
-2 (d^) (d k a) + (d q G qn ) (dLG lk ) G nk + G mn (d m a) (d n a) } 

QklQmm (^-1) (^-1) g a b + Q mk (^-l) Q nl (^-1) ^ 



tra m r m = e^Tr (fi m r m - fi m G m n9 n a) 

= e" CT Tr (-0 m (9 n G" m ) + fi m G m "9 n a) 

= (55) 

Hence we obtain 

tx£ = e-{G fc/ G m " (5 fe C) ^ " G mfe (^C) ^ (^J,) ^ ab }. (56) 

Comparing with for A 2 regarded as constant cutoff of Ag, we can write this as 
J d A y tr £ = | A (-2 - (d p 6 pa ) G* k (0,6^) + G mn (d m a) (d n a)) 

c - J d A y (-2R[G\ + G mn d m ad n a) , (57) 



assuming on-shell geometries (|18D in the last line. This formula applies for Dirac fermions, 
and with an additional factor | for Weyl fermions. It is remarkable that tr £ is essentially 
given by the appropriate curvature scalar and a contribution from the dilaton-like 

scaling factor p = e~ a . This is a very reasonable modification of the standard result fl36[) , 
as desired. 
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5 Relation with gauge theory on M.g 

We now want to interpret the fermionic action (|TT| ) as action for a Dirac fermion on the 
Moyal-Weyl quantum plane Wg coupled to a U{1) gauge field in the adjoint. This point 
of view is obtained by writing the general covariant coordinate resp. matrix Y a as 

Y a = X a + A a . (58) 

Here X a are generators of the Moyal-Weyl quantum plane, which satisfy 

[X a , X h ] = %6 ab , (59) 

where 9 ab is a constant antisymmetric tensor. These are particular solutions of the equa- 
tions of motion (JT7|). The effective geometry (|6|) for the Moyal-Weyl plane is indeed flat, 
given by 

9 ab = ac 9 bd g cd 
g ab = pg a \ det^ fe =l 
p = (det^r^H^r^c- (60) 

Consider now the change of variables 

A a (x) = -0 ab A b {x) (61) 

where A a are hermitian matrices, interpreted as smooth functions on M|. Thus we can 
write 

[Y a , f] = [X a + A\ /] = i6 ab (^Lf + i [A b , f] \ = iQ ah D h f, (62) 



giving for the quadratic form 

SVare = (2vr) 2 Tr ¥ lalb [Y% [Y b , v]/] ] 

- / d A xp^ lalb e am e bn D m D n m 



d 4 x&fl 2 A $. (63) 

This is an exact expression on Mg, where 

W 2 a = -plal h e am e bn D m D n = - TT D m D n , (64) 

and 

j a = (detg ab ) l s lb ba , {T,l b } = 2~g ab - (65) 
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We now want to rewrite the geometrical results of Section f| in terms of gauge theory on 
M.q in x-coordinates. To do this, note that most formulas of Section ^ are not generally 
covariant, but only valid in the preferred y-coordinates defined by the matrix models 
where g ab = 5 ab resp. g ab = r] ab . ( plf) defines the leading-order relation between y and x 
coordinates, 

y a = x a - e ab A b + 0(6 2 ) , (66) 



with Jacobian 



dy a 



dx b 



r)A 

i _ £~£2£ + (e 2 ) 

dx a K ' 

l-U mn F mn + 0{6 2 ). (67) 



See [4] for details of this change of variables. In order to avoid confusion we will denote 
all x-dependent tensors with a bar, and we write 

= B -^- (68) 

The Poisson tensor can be written in terms of the u(l) field strength as 

iO a \y) = [Y a , Y b ] = i6 ab - t6 ac 6 bd F cd , (69) 
where F cd is a rank two tensor in x coordinates on M.g. This amounts to 

Kb = Kb ~ F ab (70) 
at leading order. We also need the metric (||) in x-coordinates, 

Q ab = ^ac _ QaiQcjp^ (gbd _ gbegdf g ^ 

= g ac (S b c + Fj db + gJ ef F fd g db + gJ de F e} g^F gh e hb ) 

= g ac (S b c + X b ). (71) 

To compute the determinant, we use 

det (1 + Xij) = 1 + trX + - {{trXf - tr (X 2 )) + 0(X 3 ). (72) 

Then 

trX = -2F m J mn - e em g m J fn F eh g h9 F gf (73) 

and 



e a = \G ab \ 1/4 = \g ab \ 1/4 (l - ^9 mn F mn + 0(9 3 )^j . (74) 
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By a straightforward application of the relations given above one can now write down 
the second Seeley-de Witt coefficient in ^-coordinates. R[G] expressed in x-coordinates 
agrees with eq. (78) in [4], as it should be. tr£ in x-coordinates is given as 

tr£ = |# a6 | 1/4 f-^ m ^^ 



+\\9a b \ 1/A 9 mn F mn d a d a e^F pq 

= -\U\ 1/4 lg ac 9 bd F ab B 2 F cd , (75) 

where 

d 2 = B a B b g ab . (76) 

In Appendix B this expansion is given explicitly for the individual terms. We have omitted 
0(A) terms from both R[G] and tr£ , which are total derivatives and do not contribute 
to the effective action. In this way, we find for the one-loop induced action 



= J d A y (a A A + a 2 A 2 + O (logA^) 



= -4r $ - — - / d 4 y^g ac 9 bd F ab d 2 F cd A 2 . (77) 
lb n z J 2 

Finally, there is a nontrivial relation between the cutoff A of the geometrical action and 
the cutoff A of the u(l) gauge theory, which follows from the identity 

S squarc = Tr ¥ lalb [Y a , [Y\ *]] = J d 4 y^W 2 G^ = J d 4 y^^W 2 A ^- (78) 
For the Laplacians this means 

W 2 d = —W 2 a- (79) 

p 

Since we implement the cutoffs using the Schwinger parameterization ( pS|) , they are related 
as follows (cf. [4]) 

A 2 = ^-A 2 . (80) 
P 

This makes sense provided p(y)/p varies only on large scales respectively small momenta 
p A, which is our working assumption. Together with (|67|) , we obtain as a final result 
for the geometric one-loop effective action expressed in terms of gauge theory on Wg 

= " 4r * + / W)^ rP ^ )U - p) wS (81) 

where p 2 = PiPjg 1 - 1 '■ This agrees precisely with the one- loop computation in the gauge 
theory point of view obtained below. Note that the last term corresponds to tr £ in (|35"D . 
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5.1 Comparison with UV/IR mixing 



In this section, we compare the geometrical form of the one-loop effective action obtained 
in the previous section with the one-loop effective action obtained from the gauge theory 
point of view. The result is of course the same, which provides not only a nontrivial 
check for our geometrical interpretation, but also sheds new light on the conditions to 
which extent the semi-classical analysis of the previous section is valid. This generalizes 
the results of [4] to the fermionic case. We find as expected that the UV/IR mixing 
terms obtained by integrating out the fermions are given by the induced geometrical resp. 
gravitational action (|32]), in a suitable IR regime. In particular, we need an explicit, 
physical momentum cutoff A. 

Using the variables and conventions of the previous section, the action ( |TTD can be 
exactly rewritten as U(l) gauge theory on K|, which in the Euclidean case takes the form 

S[V] = (27r) 2 Tr^ t 7a [^ a ,^] 

= J ' d 4 x¥ir(B a * + ig[A a ,*}) (82) 

We introduce an explicit coupling constant g, and define a rescaled fermionic field 

%=\g ab \Te^ (83) 

in order to obtain the properly normalized effective metric g ah \ we will omit the tilde on 
\l/ henceforth. Recall also that only U(l) gauge fields are considered here, because only 
those correspond to the nontrivial geometry considered in the previous section. 

We need the 0(A 2 ) contribution to the one- loop effective action obtained by inte- 
grating out the fermionic field While this computation has been discussed several 
times in the literature [17,18,20-22], the known results are not accurate enough for our 
purpose, i.e. in the regime p 2 ,A 2 < A 2 NC where the semiclassical geometry is expected 
to make sense. We need to analyze carefully the IR regime of the well-known effective 
cutoff A e ff(p) (|8^) for non-planar graphs as p — > 0, keeping A fixed. In this regime the 
non-planar diagrams almost coincide with the planar diagrams, and the leading IR cor- 
rections due to the non-planar diagrams correspond to the induced geometrical terms in 
(32). This has not been considered in previous attempts to explain UV/IR mixing, e.g. 
in terms of exchange of closed string modes [23,24]. 

To proceed one can either square the Dirac operator as in [21], or use directly the 
fermionic Feynman rules. We choose the latter approach here, and consider the Feynman 
diagram in figure [T] corresponding to 

I\ = -^Tr\ogA -^Jd A xp^ a [A a ^}Jd A yp^ b [A b ^} 

= -iTrlogA + r*(A). (84) 
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Figure 1: fermionic one-loop diagram 



The minus sign in front is due to the fermionic loop. This gives 

n - '' d 4 k 2k a k b + k a p b + p a k b - g ab k(k + p) 



"V/ ^A a ,(p)A b ,(-p)~g a ' a ~9 b ' b 



(2tt) 4 (k ■ k)((k + p) ■ (k + p)) 
(, :! " r " - l) (85) 

which is quite close to the bosonic case, using the notation 

/c * /c = kikj cj ^ . 
k 2 = kik jg ij . (86) 

The momentum integrals are understood to be regularized by a cutoff A, implemented via 
a Schwinger parameter as in [4]. To evaluate this loop integral, we rewrite it in a different 
way as in [21] 

f d 4 k Ak a k b + 2k a p b + 2p a k b - 2g ab k(k + p) _ ik . 6 a v . _ ■ 
J (2tt) 4 ( k .k)((k+p)-(k+p)) [C ) 

d 4 k / (2k a + p a ) (2k b + p b ) - (p a p b - g ab p ■ p) 



I 



(2tt) 4 V (k- k)((k+p) ■ (k+p)) 



9ab (k-k + (k+p)-(k+ P ))Y e tkt6Jpj 1} 

f d4k ( (2k a +p a ){2k b + p b ) _ g ab \ ik . e a v . _ , 
J (2tt) 4 \(k- k){{k + p) ■ (k + p)) k-kr ' 

+ l» - *■* MttTO 11 '*"" - x) (87) 

where we replaced ^ k+p ^ k+p ^ by under the integral (which does not make a difference 
in the regularization used here). Now the first term is precisely the induced action T$ 
obtained by integrating out a scalar field $ [4], which is known to be gauge invariant. 
The second term is logarithmic and manifestly gauge- invariant. Therefore 

/d 4 
A a ,(p)A b ,(-p) g a ' a 9 b ' b (PaPb - g ab p ■ p) 
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(2tt) 4 (k ■ k)((k + p) ■ (k + p)) 

A a ,(p)A b ,(-p) g a ' a 9 b ' b (PaPb - ~9abP • p) 



for Dirac fermions, where 



1 A + iF" 



is the "effective" cutoff for non-planar graphs, and Ajvc is defined in (|60|). For the standard 
evaluation of the k— integration see e.g. [4]. To proceed we consider the IR regime 

v 2 A 2 

TT"<1. (90) 

iv JVC 

see figure 0. Then both A and A e ff are large, and we can use the asymptotic expansions 




Figure 2: relevant IR regime of A e ff(p) 




The only approximation here is the expansion (^l|) of the Bessel functions in (j88| ). T$ is 



the 1-loop effective action for a (hermitian) scalar field as computed in [4], 
g 2 1 f d 4 p / 1- . ,~ a ' a ~y h , , A 2 



2 16.* /^(-^-W^^*«" k «^ / 



2 

+i(£IF( P ))(eF(- P )) (aJ„ - • P Agjy + (47 - 30 log(^)) 



(93) 



These expressions are valid in the IR regime (EKJ) pA< A 2 NC corresponding to "mild" 
UV/IR mixing. This is the same condition which was obtained for the bosonic case [4]. 
We can then use the expansions 

ti„ = A 2 -P 2 ^ + ..., (94) 

K ff = A 4 - P 2 -^ + ..., (95) 

A 2 1 p 2 A 2 , N 

log(^) = 4^- + - (96) 

iv e// ^ NC 

which gives 

1 9 2 f d 4 p , Vi - - .p 2 A 2 



r, + 4r $ ~ \^J^9 a ' a f b F ab ( P )FM-P) 



A 4 ' 



' ,r> (/V ' W^W^^), (97) 



4 16vr 2 J (2tt 

where p 2 = p a Pb9 at '• There are obvious modifications due to the appropriate expansion of 



A 2 jj if one approaches the border of the IR regime fl9"0"|). 

To compare this with the geometrical results, we must take into account the different 



regularizations used in the heat-kernel expansion (28) and in the above one-loop com- 
putation. It was shown in [4] that these regularizations agree if we replace A 2 with 2A 2 
in the one-loop computation aboveQ. We then find complete agreement with the result 



(pl|) obtained using the geometrical point of view. Notice in particular that the induced 
gravitational action is nontrivial even in the case of e.g. N — 1 supersymmetry. This is 
now understood in terms of induced gravity, and full cancellation is obtained only in the 
case of N = 4 supersymmetry. This will be discussed below. 

Finally, and T$ can be related directly to the geometrical induced action (|32|) in a 
more restricted IR regime, as in [4]. Assume first that 



6 while this was strictly speaking established only for the bosonic case, the argument should extend to 
the fcrmionic case without difficulties. 



19 



A < A NC - 

Then the IR regime (|90| ) amounts to 

P < ^nc, (98) 
which is very reasonable range of validity for the classical gravity action. In this case, 



A^ = -^AV « AY ~ K'p.p (99) 



so that we can replace 



A e// - ^ " P A lff ~ A 4 - ip • P A 2 . (100) 
Then the leading contribution to T$ is 

r. ~ 4 / 1^ (t™^^^ 

-^(0F(p))(6F{-p))p.p +0(log(A)) +finite) 
+0(log(A)) + finite), (101) 



where again p 2 = p a Pb9 ab - Taking into account again the appropriate replacement A 2 — > 
2A 2 corresponding to the geometrical regularization in (p8~ ) and setting = 2 for Dirac 



fermions, one finds as in [4] complete agreement between the above result for with the 
result (|32|) obtained from the geometrical point of view. 

Assume finally that the condition A Ajvc is violated, while maintaining the IR 
regime (|90|). Then there are additional terms in the effective action T$ beyond ( |101| ). 



They correspond to noncommutative corrections beyond the semi-classical geometrical 
terms in (P^D; we refer to [4] for some explicit results. 



Gauge fields resp. gravitons The complete one-loop effective action for the gravita- 
tional sector sector of the basic matrix model ([!]) contains contributions of scalar fields, 
fermions, and "gauge fields". The latter means trace- U(l) components resp. gravitons 
inherent in G a b, as well as possible nonabelian SU(n) components [1]. Ignoring nonabelian 
fields for now, it remains to compute the contributions of the U(l) gauge fields in a loop 
with 2 external A- legs. This would lead to a similar contribution, denoted by 

e~ TA = / dAe~ s . (102) 

J one— loop 
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While Ta is essentially straightforward to compute following e.g. [21], we can short-cut 
this computation by taking advantage of the finiteness of the N = 4 supersymmetric 
extension of the model. Note that the presence of particular scalar interaction terms in 
the N = 4 model is irrelevant for ( |102| ) at one loop. This leads to the result ( |103| ) given 
below. 



5.2 Cancellations and super symmetry 

It is very interesting to compare the fermionic and the bosonic contribution to the gravi- 
tational action. As is well-known [17,21], we note that the fermionic contribution to the 
one-loop effective action in NC gauge theory does not quite cancel the scalar contribu- 
tion, due to (|92|). This means that even in supersymmetric cases some UV/IR mixing 
may remain. From the geometrical point of view, this terms corresponds to a gravita- 
tional action tr£A 2 = — 2R[G] A 2 + so that the cutoff A 2 should be interpreted as 
effective gravitational constant This is completely analogous to the commutative case, 



where the gravitational term ( |5B| ) is induced. The remaining UV/IR mixing term cancels 
only in the case of N = 4 supersymmetry. We can therefore identify A as the scale of 
N = 4 SUSY breaking (assuming such a model), above which the model is finite. These 
observations strongly suggest that for the model to be well-defined at the quantum level, 
N = 4 SUSY is required above the gravity scale i.e. the Planck scale. This is realized by 
the IKKT model [5] on a NC background. 

Let us therefore consider an extension of the matrix model (P with ng scalar fields 
(= hermitian matrices (p l ) and Dirac fermions (hence n$ = ~ for Weyl fermions). The 
model with N = 4 SUSY has n$ = 6 and = 2, in addition to the U(l) gauge field 
resp. the graviton h a b- Taking its finiteness for granted, it follows that 

r A = -2r*-6r$. (103) 

This also holds for the SU(n) components of the gauge fields in a nonabelian versions of 
this N = 4 model; cf. [1]. Using fl35|) , this can be written explicitly as 

Ta = t^J dAy (~ 4 ^ 4 " { l R[d] + 2tl£) ^ 2 + °( logX )) • ( 104 ) 

This almost coincides with the contribution of 2 scalars fl34]), except for tr£ which stands 

for (0). At this point, it is important to keep in mind that a term J d 4 y \J \G\ A 4 in 
the framework of emergent gravity does not amount to a large cosmological constant, as 
discussed in [1,4]; see also [25]. We only point out here that flat space is a solution even 
in the presence of this term. These issues will be discussed elsewhere in more detail. 



N = 4 supersymmetry transformations. The explicit form of the A^ = 4 SUSY 
transformations can be obtained as follows. Consider the basic matrix model (|l|) in 10 
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dimensions with fermions, with action 



S = -Tr \^-[Y a , Y b ] [Y a \ Y b ')r) aa ,r) W + * la [Y a , *] J , a = 0, 9, (105) 

where \I/ is a 10- dimensional Majorana-Weyl spinor. It is well known that this action 
enjoys the 2 "matrix" supersymmetries [5] 

§ (i)q = l -[Y a ,Y b ] lab e, 5^Y a = ie 7 a ^, 

<5 (2) ^ = f, 5 {2) Y a = (106) 

where = \ Yf a , 7b], and e, £ are Grassmann- valued spinors. To recover spacetime super- 
symmetry, we split the matrices into 4 + 6 dimensions as 

Y a = (Y^^), /i = 0,...,3, i = l,..,6. (107) 

Then the 4-dimensional Moyal-Weyl quantum plane Mg is a (BPS) solution of the gener- 
alized matrix equations of motion, embedded as 

Y" = X", fi = 0,...,3, 

ft = 0. (108) 

All previous geometrical considerations can be generalized, except that the matrix model 
now contains scalar fields <f) l (y)- Even though we could consider a general metric G a b as 
above, let us for simplicity only discuss the point of view of U(l) gauge theory on M.g. If 
we set £ = |6> Ml/ 7^e following [5] and use 7 a (3+;) = 7a73+i recalling Y 11 = X^ — 9^ h ' A v and 
QS5D) then the combined transformation 5 = 5^ + 6^ takes the form 

5A U = -iVpguvfiT* (109) 

where £^ = | [7 / \7 1/ |- The constant factors of p can be absorbed by rescaling the fields. 
Noting that 73+4 = 75 Aj where Aj define the 6-dimensional Euclidean Clifford algebra, 
this indeed amounts to (abelian) N = A supersymmetry on R|. 

In the case of general NC backgrounds, the SUSY transformation will also act on 
the metric G liV {y). This can be viewed as a supersymmetric form of emergent gravity, 
which will be worked out elsewhere. It also means that nontrivial background geometries 
"spontaneously" break N = 4 SUSY, as desired. 
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6 Discussion and outlook 



In this paper, fermions are studied in the framework of emergent noncommutative grav- 
ity, as realized through matrix models of Yang-Mills type. The matrix model strongly 
suggests a particular fermionic term in the action, corresponding to a specific coupling to 
a background geometry with nontrivial metric G^. This coupling is similar to the stan- 
dard coupling of fermions to a gravitational background, except that the spin connection 
vanishes in the preferred coordinates associated with the matrix model. 

The main result of this paper is that in spite of this unusual feature, the resulting 
fermionic action is very reasonable, and properly describes fermions coupled to emergent 
gravity. In the point particle limit, fermions propagate along the appropriate trajectories, 
albeit with a different rotation of the spin. At the quantum level, we find an induced 
gravitational action which includes the expected Einstein-Hilbert term with a modified 
coefficient, as well as an additional term for a scalar density reminiscent of a dilaton. There 
are further terms which vanish for on-shell geometries. We conclude that the framework of 
emergent gravity does extend to fermions in a reasonable manner, and might well provide 
- in a suitable extension - a physically viable theory of gravity. 

In a second part of the paper, we compare this induced gravitational action with the 
well-known UV/IR mixing in NC gauge theory due to fermions. Generalizing the results 
in [4] for scalar fields, we find as expected that the UV/IR mixing can be explained pre- 
cisely by the gravitational point of view. This also provides a nice understanding for 
the fact that some UV/IR mixing remains in supersymmetric cases, and only disappears 
for N = 4 supersymmetry. The reason is that a gravitational action is induced even in 
supersymmetric cases, except in N = 4 SUSY. This in turn leads to the conjecture that 
the gravitational constant should be related to the scale of N = 4 SUSY breaking, which 
is quite reasonable. All of these findings suggest that the IKKT model on a noncommu- 
tative background [5-8] should be the most promising candidate for a realistic version of 
emergent gravity. These issues will be discussed in more detail elsewhere. 
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7 Appendix A: Computation of R 

We quote some identities which appear in the computation of R[G] and R[G] in terms of 
6*-vielbeins: 

(d b G bd )G ca (d a G cd ) = -G ac {d c e^)G bd {d d e- q 1 )g^-2{d a e^)G bc {d c e^) 
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G ab G dbddGac = 

(d b G bd )(d d a) = 

G™(d p G mn )(d q G mn ) = 

G np (d p G ac )(d c G na ) = 



- (d a 9 ap ) (d b 9 b «) g pq (110) 
2G mp G n W m 1 a d p d q 6- b 1 g ab + G bc {d^) G ad {d d 6 bq l ) g™ 

+G ad (0,9-;) G bc (d c 9^) <p (111) 

-2rG"WJ, + 2G mn G pq (d p 9-l) (d q 9j) g ab (112) 

G mn (d m a) (d n a) + 9 bc (d b 9 da ) (d d a) g ca (113) 
-2G« (d p 9 mn ) (d q 9^ n ) 

-2G mn G pq (d p 9 ma ) (d q 9~l) g ab (114) 

_ 2G m P {dp9 n q) (gj-lj _ {dJCq) {dc9 n P) ^ 

-G nk (d.9- 1 ) G d (d^) g 



G mn d m ad n a 



G 



mn d d a 



_ GPqg mn (^-1 j qM fag-la) 

\g* (d p 9™) [0,9-1) + \e mn G pq d p d q ^ n 



(115) 



Below are some identities which have not appeared in the Ricci scalar but appear in the 
computation of tr£: 



G mn (d m G np ) g ab 9 qa (d q 9 pb ) 



d m d n G n 



G nk 9 qa [d q 9 nb ) 9 lc {d t 9 kd ) g ab g cd 
G nk 9 qa (d q 9 nb ) 9 lc [d l 9 kd ) g ac g bd 
G nk 9 qa (d q 9 nb ) 9 lc (d^) g ad g bc 
G mn (d p G pm ) (d q G qn ) 



-G mk (d k 9- a ) G nl [dtf-l) g ab - (d m 9 ma ) G nk {d k 9~l) 
29 ma d m d n 9 nb g ab + (d m 9 na ) (d n 9 mb ) g ab 
+ (d m 9 ma ) [d n 9 nb ) g ab 
G mk (d k 9 ma ) G nl {d l 9~J) g ab 

QklQmn (^-l) (^-l) g ab 

G mp [d p 9~l) G nq (d q 9 m l) g ab 

(d p 9 pa ) (d q 9 qb ) g ab + 2 (d p 9 pa ) G qk (d k 9^) 

+G mk (d k 9 ma ) G nl {di9~l) g ab (116) 



8 Appendix B: Expressing R and trS in x coordinates 

Let us rewrite the terms which compose R and tr£ in terms of the u(l) gauge fields. We 
will need [4] 

e° = (det G ab ) 1/4 = (det g ab ) 1/4 (l - l -9 mn F mn + 0{9 2 )^ , 

G ab = ^ac _ (jaiQcjp^ (gbd _ gbegdf j?^ g ^ 
gab gab Qacgbdp 

dab' = Kb 'Fab- (H7) 
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and denote | det g a b\ = \g a b\- For the following terms this gives to 0(A 2 ) in x-coordinates 



d'ye-* {d m 6 ma )G nk (d k 9^) 
d A ye-° (d m 6 na ) (d n 6 mb ) g ab 
d A ye-°9 ma d m d n 6 nb g ab 



d 4 y e _cr G pq (d p 6 mn ) [d q 9^ n ) 



Ae- ff G H G™ {d k e- m l) {dtf-Dg^ 



d%e-°G mk {d k 6-l)G nl {d l 9^)g ah 
where we used [4] 

J d 4 xF m J mn d a d a F pq 6™ = aJ d 4 xg fh F fa d a d s F hs , 
d 4 xg ab F am d m d s F bs = 



-\ J d 4 x\-g ab \ l i A 9 mn F mn d a d a ep% 

-i J d 4 x\g ab \ l l A e mn F m J) a d a e^F % 
- l - J d A x\g^e mn F mn &dJ p % 
- j d'xlg^^BJ^F^ 



d 4 x\g ab \ x ' 4 d mn F mn d a d a P% 



VI 



pq 



+-e mn F mn 8 a d a e™F m 

j d 4 x\g ab \ 1 ^(-~6 mn F mn d a 8J™F pq 

-\g ac g hd F ab d 2 F cd ^ 

-\ [ d A x\g ab \ l ^6 mn F mn d a d a 6™F pq (U8) 



1 



d 4 x ^g ah 9 mr F hm B 2 F ra - g si g ah F hs d l diF ia . (119) 



9 Appendix C: a k and the spin connection. 

We observe that using the Jacobi equation, the linear term a k appearing in I/) 2 
be written in the form 

a k = G kl d l a+- A [ la , lb ]6 lk {d l 6 ab )-T k 



[7a, 7fc] lk {d^) + 6 ka G mn (d n 9 m l a ) 
haM lk {di6 ab ), 



can 



(120) 
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where the last line holds only for on-shell geometries, using the equations of motion. In 
the standard form of the Dirac operator, this linear term would have the form 



[ 7a , 7k] ^i b G lk - r 



where 



UJ 



a b 



7 [7a, 76] kd g dm {d c 6 mb ) g™ + -6 lk {d^) 



-9r c l g ca V k 9 lb 

-i„ma alb 



er c L 9 ca d k e lb - e- c l g ™T km e mh 

6 lb {d^ 1 ) <T - 9 na (d k 9~l) g 



l )g cb 

1\ „cb 



-e mb {d m e^) <r + e ma {d m e~ k l c ) g 

nmb ( o ana\ r~i , nma I a Qnb\ r~i 



G kl d l a, 



-UJ 



b a 



:i2v 



:i22) 



is the spin connection, using the explicit form of the frame (0). Note that the effect of 
the spin connection is different in our case. Nevertheless, it turns out that tr£ can be 
rewritten as Ricci scalar R[G] plus a term dependent on a, as shown in Section P~2. 
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